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摘 ， 要 :工程 尺度 的 溶 蚀 岩 体 难 以 开展 室内 外 力学 试验 ,导致 溶 蚀 岩 体 的 力学 参数 获取 存在 困难 ， 
因此 基于 岩石 矿物 含量 特征 、 室 内 力学 试验 和 岩 体 结构 面 特 征 ,采用 三 维 颗粒 流离 散 元 法 ,通过 建 
立 等 效 孔 阶 型 溶 蚀 岩 体 模 型 进行 单 轴 压缩 数值 试验 ,进而 分 析 其 力学 特性 和 变形 破坏 机 制 。 研 究 
表明 :采用 改进 平行 黏 结 模型 ,基于 细 观 矿物 特征 ,通过 元 胞 自动 机 算法 吻 除 颗粒 可 建立 孔隙 型 溶 
蚀 岩 体 模 型 ;加 载 前 期 , 岩 体 结构 面 首先 快速 破坏 而 产生 以 前 切 为 主 的 微 裂纹 ; 随 加 载 进行 , 岩 块 内 
逐渐 产生 以 拉 裂 纹 为 主 的 破坏 ,其 微 裂 纹 呈 指数 增加 ,而 结构 面 微 裂纹 先 激增 后 趋 于 稳定 ;相同 轴 
向 应 变 时 , 岩 块 内 拉 裂 纹 随 溶 蚀 率 增 加 而 增加 ,而 结构 面 剪 切 微 裂 纹 减 少 ; 岩 体 的 变形 破坏 分 为 结 
构 面 快速 破坏 、 岩 块 弹性 变形 、 岩 块 塑 性 变形 和 岩 体 完全 破坏 等 4 个 阶段 ,其 破坏 形式 随 溶 蚀 率 增 
加 而 从 整体 均匀 性 破坏 转化 为 局 部 结构 性 破坏 ; 溶 蚀 使 岩 体 强度 降低 , 溶 蚀 率 与 单 轴 抗 压强 度 和 变 
形 模 量 分 别 呈 反比 和 负 指 数 函 数 关系 。 

关键 词 : 溶 蚀 岩 体 ; 离散 元 法 ; 单 轴 压缩 ;力学 特性 ; 岩 体 强度 

中 图 分 类 号 :TU452 ”文献 标志 码 :A DOI:10.11776/j. issn. 1000-4939. 2023. 05. 020 


Numerical test on the mechanical properties of dissolution rock 
mass under the uniaxial compression condition 
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Abstract: lt is difficult to carry out laboratory and field mechanical tests on engineering-scale dissolution 
rock mass , which leads to great difficulty in obtaining mechanical parameters of the dissolution rock mass. 
Therefore ,based on the characteristics of rock mineral content, laboratory mechanical tests ,and rock mass 
structural plane characteristics ,three-dimensional particle flow of discrete element method was adopted to 
conduct a uniaxial compression numerical test by establishing an equivalent dissolution rock mass model. 
Finally ,the mechanical properties and deformation failure mechanism are analyzed. Based on the charac- 


teristics of meso-minerals and macro-structural planes ,the research shows that the establishment of an e- 
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quivalent dissolution rock mass model can truly reflect the mechanical behavior of it. In the early stage of 


loading ,the structural plane of rock mass was quickly destroyed to produce micro-cracks mainly of shear. 


However ,with the loading progressing ,tensjle failure occurred gradually in the rock mass ,and its micro- 


cracks increased exponentially ,while the micro-cracks of the structural plane surged first and then stabi- 


lized. Under the same axial strain ,the internal tensile cracks in the rock block increase by the dissolution 


rate ,but the shear microcracks on the structural plane decrease. The deformation and failure of the rock 


mass can be divided into four stages :rapid failure of the structural plane ,elastic deformation of the rock, 


plastic deformation of the rock ,and complete failure of the rock mass. The failure form of the rock mass 


transitions from the overall-uniform failure to the local-structural failure with the increase of the dissolution 


rate. dissolution reduces the strength of rock mass. The dissolution rate is inversely proportional to the uni- 


axial compressive strength ,and has a negative exponential relationship with the deformation modulus. 


Key words:dissolution rock mass; discrete element method ; uniaxial compression ; mechanical properties ; 


strength of rock mass 


岩 体内 节理 、 裂 隙 等 地 质 缺 陷 对 其 力学 响应 特 
性 一 直 是 工程 界 关注 的 重点 和 难点 '“  。 结 构 面 对 
岩 体 力学 性 质 起 控制 作用 中 ,而 在 约 占 全 球 地 表面 
积 1/5 ,全 国 面积 17/3 的 喀斯特 地 区 ,一 定 深度 范围 
内 发 育 着 不 同 深 蚀 程度 的 岩 体 “ ,其 力学 性 质 不 仅 
受 结 构 面 控制 ,还 受到 大 量 溶 蚀 和 孔 际 、 裂 际 的 影响 ， 
使 得 溶 刨 告 体 的 力学 特性 更 加 复杂 。 目 前 ,关于 市 
理 . 裂 际 岩 体 的 力学 特性 研究 取得 了 丰富 的 成 
果 " ,而 关于 溶 蚀 岩 体 力学 特性 研究 不 多 。 随 
着 工程 建设 不 断 发 展 , 溶 蚀 岩 体 地 基 和 边 坡 稳定 性 
问题 日 益 显 著 “ 2 ,喀斯特 地 区 的 岩溶 缺陷 对 岩 体 
力学 特性 的 影响 吸 待 解决 。 现 阶段 , 溶 蚀 岩 体 研究 往 
往 只 关注 个 别 较 大 岩溶 洞 隙 的 影响 ,而 忽视 数量 庞大 
的 中 小 尺度 告 溶洞 阶 对 岩 体 力学 特性 的 作用 。 

近年 来 ,数值 模拟 分 析 成 为 研究 溶 蚀 宕 体力 学 


符合 实际 溶 蚀 发 育 特 征 的 岩 体 模型 ,进而 对 溶 蚀 岩 
体力 学 特性 进行 研究 。 

工程 尺度 的 岩 体 现场 力学 试验 实现 较 困难 且 费 
用 高 ,现场 岩 体 力学 原型 试验 难以 普及 ,而 数值 分 析 
方法 对 溶 蚀 岩 体力 学 特性 研究 逐渐 被 认可 。 离 散 元 
颗粒 流 法 能 允 近 岩 土 材 料 的 力学 相应 特性 ,在 岩 体 
力学 特性 方面 的 研究 方面 具有 很 强 的 优势 "1。 部 
分 学 者 在 采用 颗粒 流离 散 元 法 对 洞 隙 影响 岩石 
( 体 ) 力 学 特性 研究 方面 取得 了 重要 进展 。 如 黎 崇 
金 等 利用 PFC 研究 大 理 岩 的 裂纹 扩展 及 力学 特 
性 ; 杨 圣 奇 等 站 利用 PFC 研究 单 也 和 双 孔 洞 阶 对 
砂岩 的 力学 特性 影响 ; 张 社 荣 等 1 采用 PFC” 从 细 观 
力学 角度 研究 不 同 的 溶 蚀 孔隙 结构 对 灰 岩 力学 特性 
的 影响 。 学 者 们 基于 二 维 离散 元 颗粒 流 , 从 细 观 力学 
角度 分 析 了 多 种 孔洞 对 岩石 力学 行为 特性 的 影响 。 


特性 的 重要 手段 ,而 研究 溶 蚀 岩 体力 学 特性 的 前 提 
是 构建 溶 蚀 岩 体 模型 。 目 前 ,在 构建 溶 蚀 宕 体 模型 
方面 的 研究 取得 了 一 定 的 成 果 , 如 张 菊 明 等 ”通过 
采用 椭 球 洞 体 建 立 随 机 洞 际 溶 蚀 岩 体 模型 ; 陈 祥 军 
等 ”针对 小 尺度 的 岩溶 洞 际 建立 溶 蚀 岩 体 随 机 结 
构 模 型 ; 刘 宏 等 “采用 坚 向 椭圆 形 溶 孔 建立 溶 蚀 岩 
体 模 型 ,进而 研究 了 溶 蚀 岩 体 地 基 稳 定性 等 力学 特 
性 ; 张 社 荣 等 ”建立 了 溶 蚀 岩 体 参数 完全 随机 场 
力学 模型 ,并 利用 Bayes 法 约束 该 模型 的 空间 变异 
性 ,进而 研究 溶 蚀 宕 体 坝 基 稳定 性 等 力学 特性 。 针 
对 中 小 尺度 溶 蚀 洞 际 , 基 于 现场 溶 蚀 洞 隙 的 统计 构 
建 随机 溶 蚀 岩 体 模型 的 方法 ,已 成 为 研究 溶 蚀 岩 体 
力学 特性 的 重要 手段 。 但 建立 的 溶 蚀 岩 体 模型 较为 
简单 ,与 实际 溶 蚀 特征 匹配 度 不 够 ,因此 有 必要 构建 


而 三 维 空间 内 ,基于 细 观 矿物 组 分 和 宏观 结构 面 特征 
的 溶 蚀 岩 体 的 细 观 力学 机 制 ,有 待 展开 进一步 研究 。 

因此 ,本 研究 以 贵州 典型 喀斯特 地 区 溶 蚀 岩 体 
为 研究 对 象 ,基于 野外 调查 、 偏 光 显微镜 和 室内 力学 
试验 等 研究 手段 ,采用 三 维 颗粒 流离 散 元 法 建立 了 
考虑 尺寸 效应 的 溶 蚀 宕 体 模型 ,进而 探讨 单 轴 压 缩 
条 件 下 不 同 溶 蚀 程度 孔 际 宕 体力 学 特性 响应 。 


1 颗粒 流离 散 元 法 


颗粒 流离 散 元 法 最 早 由 CUNDALL 等 25 提出 ， 
该 方法 将 岩 体 离散 成 刚性 细 观 颗粒 (但 力学 关系 上 
人 允许 重 麦 ) ,基于 牛顿 第 二 定律 定义 颗粒 间 的 力学 运 
动 行为 ,采用 接触 (contact) 模型 来 反映 颗粒 间 的 力 
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与 位 移 关系 。 目 前 ,在 岩 体力 学 研究 领域 ,线性 接触 
模型 .平行 竺 结 模型 和 光滑 节理 模型 是 广泛 运用 的 
3 种 接触 模型 。 


1.1 改进 的 平行 猪 结 模型 


平行 黏 结 模型 ( parallel bonding model, PBM ) 是 
通过 平行 条 结 将 相互 接触 的 颗粒 黏 结 成 集合 体 以 表 
征 完整 的 岩 块 ,可 看 作 一 组 具有 恒定 法 向 和 切 向 刚 
度 的 弹簧 ,均匀 地 分 布 在 颗粒 间接 触 面 上 。PBM 能 
够 承受 拉 应 力 且 能 抵抗 颗粒 相对 旋转 ,平行 黏 结 
坏 时 ,退化 为 线性 接触 模型 ”” 。 

平行 慕 结 力 F 与 力矩 M 分 成 法 向 与 切 向 分 


量 , 即 


F=Fn+rFt 
| (1) 


M= Mn+Mt 
式 中 :FF JW .分别 为 力 与 力矩 的 法 向 与 切 向 
分 量 大 小 ;nt 分 别 为 法 向 与 切 向 单位 向 量 。 在 三 维 
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拉 问 题 。 定 义 B 为 软化 刚度 天 与 荷载 刚度 无 之 比 ， 
即 B = /i ,其 中 = 和 4 。 其 中 , 当 B 一 w 时 , 即 
为 传统 完全 脆性 的 PBM 模型 (图 1) 。 


A 
图 1 软化 平行 竺 结 模型 的 法 向 力 -位 移 法 则 
Fig. 1 Normal force-displacement law of softening PBM 
平行 黏 结 的 拉 伸 法 向 力 F, 达 到 并 超过 法 向 笑 
结 强度 极限 时 发 生 软 化 , 即 软化 临界 为 
Fw =— 0.4 (6) 


软化 后 , 颖 结 的 法 向 力 与 位 移 仍 为 线性 关系 。 


颗粒 流 中 , 切 向 分 量 法 向 力 F 和 力矩 村 , 沿 接 触 平面 
5-t 坐标 系 分 为 FF, 和 MM.。 

颗粒 间 的 相对 位 移 及 转动 会 引起 弹性 力 与 力矩 
增 量 ,每 一 个 计算 步 中 , 力 与 力矩 增 量 》 


一 = -有 4A56 ,AP = -大 4A6. 


ER 2 (2) 
AM, = -JIJAO ,AM = -大 17AO、 

式 中 :A5, .A6, 和 Ab,、Ab. 分 别 为 法 向 、 切 向 位 移 增 

量 和 相对 扭矩 . 弯 矩 增 量 , 其 中 ,法 向 力 拉 为 负 , 即 法 

向 位 移 6, >0 为 拉 伸 ; h,、k. 分 别 为 平行 黏 结 的 法 

向 与 切 向 刚度 ,由 式 (3) 求 出 ; 4 、 六 了 分 别 为 平行 

黏 结 接触 面积 .惯性 矩 和 极 惯性 矩 , 由 式 (4) 求 得 。 
k, = E/(R, +RR) 有 = 大] (3) 


4 = nh’, 1 = mR /4,] = 7R'/2 (4) 
式 中 : 及 =Amin(R ,Rn ) ,R。 和 Rs 为 两 个 相 邻 接触 
颗粒 半径 , 和 为 半径 乘 数 。 
平行 慕 结 接触 面 的 正 应 力 o 和 剪 应 力 7 为 


二 三 |- 
和 = 一 +- 一 
4 了 
(5) 
RD 
4 4 


传统 的 PBM 模型 很 难 同 时 满足 实际 岩 体 较 高 
的 抗 压 与 抗 拉 强度 之 比 ,而 文献 [29 ] 在 平行 苛 结 模 
型 的 法 向 分 量 中 引入 软化 系数 B, 能 很 好 解决 高 压 


由 于 法 向 分 量 发 生 软 化 ,破坏 判 据 由 法 向 恭 结 强度 
变 为 法 向 位 移 , 即 满足 式 (7) 时 发 生 拉 伸 破 坏 。 


6, + RI 9.1 > Su (7) 
式 中 ,5 是 软化 黏 结 最 大 位 移 量 ( 拉 伸 ) , 即 
_ Cef1+ 
Or | (8 ) 


而 切 向 分 量 不 发 生 软化 ,其 剪 切 黏 结 强度 r. 遵 
循 库仑 准则 。 当 7 > r。 时 , 黏 结 被 剪断 ,产生 剪 切 列 
纹 。 当 平行 儿 结 发 生 剪 切 或 拉 断 破坏 时 ,模型 退化 
为 只 有 线 弹 性 和 摩擦 作用 。 


1.2 光滑 节理 模型 


光滑 节理 模型 (smooth joint model,SJM ) 通过 模 
拟 节 理 两 侧 一 定 范围 内 的 接触 特性 来 反映 前 胀 . 摩 
擦 和 条 结 节理 的 力学 特性 。 节 理 接触 部 位 是 光滑 
的 ,两 侧 颗粒 可 以 相互 覆盖 或 滑动 ,不 同 于 黏 结 模型 
(PBM) 中 一 个 颗粒 绕 过 另 一 个 颗粒 。 当 节理 秋 结 此 
裂 时 ,光滑 节理 模型 退化 为 只 具有 线 弹 性 和 摩擦 
特性 。 

在 创建 光滑 节理 模型 后 , 黏 结 模型 转 为 光滑 节 
理 模 型 ,光滑 节理 模型 中 多 数 参 数 继承 笑 结 模型 的 
参数 。 光 滑 节理 模型 可 看 作为 一 系列 弹 得 均匀 分 布 
在 一 个 颗粒 间接 触 面 ”。 

在 每 一 个 计算 步 中 ,将 位 移 增 量 分 量 分 别 与 节 
理 接触 面 法 向 与 切 向 刚度 相 乘 即 可 算得 各 力 的 分 量 
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增 量 , 则 力 分 量 更 新 为 
F = 下 + 天 4A6， 
| =F -EAAS. 
节理 的 切 向 黏 结 强度 ~ 遵从 库仑 准则 ,由 黏 聚 
力 cc 和 内 摩擦 角 四, 及 正 应 力 r(c = 了 /4) 决 定 。 
有 具 黏 结 特性 的 节理 ,接触 面 法 向 力 小 于 节理 法 
向 黏 结 强 度 极 限时 ( 即 ,< -4) ,节理 黏 结 发 生 
拉 伸 破坏 ,产生 拉 伸 微 裂 纹 ; 而 接触 面 切 向 应 力 超 过 
抗 剪 强度 极限 时 ( 即 ,> r.4) , 黏 绪 发 生 剪 切 破坏 ， 
产生 剪 切 微 裂 纹 。 


(9) 


2 颗粒 流 细 观 参数 标定 


为 了 颗粒 流 试验 能 够 模拟 岩 体 的 力学 行为 , 需 
要 通过 不 断 改变 细 观 参数 进行 与 室内 力学 试验 相 匹 
配 的 颗粒 流 数 值 试 验 , 即 细 观 参数 标定 。 


2.1 灰 岩 细 观 结构 特征 


由 于 矿物 晶体 之 间 的 胶结 强度 各 不 相同 ,采用 
统一 的 细 观 参数 构建 所 有 颗粒 间接 触 模型 ,难以 反 
映 不 同 矿物 颗粒 组 合影 响 岩 石 强度 的 差异 。 因 此 ， 
非常 有 必要 基于 细 观 矿物 组 合 特征 , 按 一 定 比 例 建 
立 不 同 颗 粒 间 的 接触 模型 。 

本 研究 选取 新 鲜 的 灰 岩 试 样 ,将 其 切割 并 打磨 
成 厚度 为 30 hm 的 薄片 ,利用 加 拿 大 树胶 将 其 平整 
固定 于 玻 片 上 而 制 成 岩石 薄片 ,再 利用 蕴 素 红 -S 染 
色 剂 对 薄片 进行 染色 (方解石 被 染 为 深 红色 ) , 置 于 
偏光 显微镜 透射 光 下 进行 鉴定 分 析 。 鉴 定 结 果 表 
明 ,岩石 试 样 为 细 唱 含 石英 灰 岩 , 呈 细 粒 结构 ,颗粒 
间 紧 密 接触 ,为 颗粒 支撑 类 型 和 以 亮 品 方解石 胶结 
物 为 主 的 接触 式 胶结 (图 2) 。 


图 2 偏光 显微镜 下 灰 岩 细 观 矿物 特征 


Fig.2 Microscopic characteristics of limestone 


under polarizing microscope 


灰 宕 主要 矿物 成 分 为 方解石 ,含量 88% ;还 有 
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少量 石英 和 黏土 矿物 ( 陆 源 物 质 ) 等 ,含量 分 别 为 
9% 和 3% ,零星 分 布 于 方解石 中 。 


2.2 室内 岩石 力学 试验 


为 确保 颗粒 流离 散 元 试验 能 准确 反映 洲 蚀 岩 体 
的 力学 行为 ,需要 获取 岩 体 基 本 力学 参数 ,以 此 作为 
细 观 标定 的 参考 依据 。 首 先 ,选取 完整 性 较 好 的 灰 
岩 制 成 50 mm x 100 mm( 直径 x 高 ) 的 标准 圆柱 试 
样 ,采用 应 力 控制 (速率 为 0.5 ~ 1.0 MPavs ) 的 方式 
在 伺服 压力 机 进行 单 轴 压 缩 试验 ,获取 单 轴 抗 压强 
度 oo、 弹性 模 量 £ 和 泊 松 比 v; 其 次 , 选 完整 性 较 好 
的 灰 岩 岩 块 制 成 30 mm x25 mm( 直径 x 高 ) 标 准 圆 
柱 试 样 ,采用 速率 为 0. 1 MPavs 的 应 力 控制 在 伺服 
压力 机 进行 巴西 臂 裂 试验 获取 单 轴 抗 拉 强 度 c,, 如 
图 3(a) 所 示 。 


轴 向 应 变 e 
0 0.001 0.002 0.003 0.004 


单 轴 压 缩 试验 


巴西 驶 裂 试 验 


E=21.55 GPa 


轴 向 应 力 o,/MPa 
拉 应 力 o,/MPa 


| 
a | 
0=11.30 MPa \# 


0 0.005 0.010 0.015 
轴 向 应 变 e 
(a) 岩石 试验 


ID 
a 


T=0.45+0.810 


最 大 剪 应 力 t,,/MPa 
f5 
入 


正 应 力 o/MPa 


(b) 结构 面 剪 切 试验 
图 3 室内 岩石 和 结构 面 力 学 试验 
Fig.3 Indoor mechanical tests of rock and discontinuity 
最 后 ,采集 含 结构 面 的 灰 岩 试 样 ,将 其 切割 打磨 
成 长 宽 高 为 50 mm x50 mm x50 mm 的 方 样 ,利用 便 
携 式 直 剪 仪 进行 剪 切 试验 ,采用 正 应 力 恒 定 , 逐 级 施 
加 剪 应 力 ,每 30 s 施加 一 级 ,直至 试 样 破坏 。 通 过 拟 
合 抗 剪 强度 - 正 应 力 关 系 曲线 求 取 结构 面 抗 剪 强度 
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参数 黏 聚 力 C 和 内 摩擦 角 p, 如 图 3(b) 所 示 。 室 内 
试验 结果 为 : 单 轴 抗 压强 度 为 171. 10 kPa .弹性 模 量 
为 21. 55 GPa 泊 松 比 为 0. 19 和 抗 拉 强度 为 11. 30 
kPa ,结构 面 秋 聚 力 为 0.45 MPa, 内 摩擦 角 约 为 39"。 
单 轴 压 缩 试验 试 件 的 破坏 形态 呈 典 型 柱状 臂 裂 脆性 
拉 破 坏 形 式 , 而 含 结构 面 直 剪 试验 岩石 试 件 呈 剪 切 
破坏 形式 。 


2.3 ”颗粒 流 匹 配 试验 细 观 标定 


2.3.1 模型 建立 
为 准确 获取 颗粒 流 的 细 观 力学 参数 ,需要 利用 
颗粒 流 匹 配 试验 进行 细 观 参数 标定 。 根 据 灰 岩 细 观 
矿物 结构 和 矿物 含量 特征 (图 2) ,将 矿物 粒 径 放 大 
一 定 比 例 后 ,采用 自动 元 胞 机 方法 ”建立 与 室内 力 
学 试验 试 样 尺 寸 一 致 的 匹配 试验 颗粒 模型 (图 4a)， 
再 将 颗粒 及 颗粒 间 赋 予 接触 模型 , 即 全 部 非 裂 际 面 
颗粒 赋予 平行 黏 结 模型 (PBM) ,而 在 用 于 剪 切 试验 
的 试 样 中 预制 一 条 贯穿 的 裂 际 面 , 裂 隙 面 两 侧 颗粒 
间接 触 采 用 光滑 节理 模型 (SIM) ,从 而 建立 了 单 轴 
压缩 匹配 试验 .巴西 辟 裂 试验 和 结构 面 剪 切 试验 的 
接触 模型 (图 4b)。 
CN 3 @ 方解石 顺和 
图 厂 英 颗粒 
全 条 十 矿物 颖 粒 


巴西 辟 裂 试验 结构 面 剪 切 试验 


单 轴 压 缩 试验 


(b) 接触 模型 
图 4 数值 匹配 试验 颗粒 及 其 接触 模型 
Fig.4 Particle and contact models of 


numerical matching test 
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2.3.2 模拟 试验 细 观 标定 
在 单 轴 压 缩 数值 试 验 中 ,采用 模型 顶 底 部 的 墙 
体 (wall) 作为 试 件 的 加 载 板 , 墙 体 与 颗粒 之 间 为 线 
性 接触 模型 ,摩擦 系数 为 0。 试验 时 ,采用 应 变 率 为 
5 x10-《s 进行 控制 ,并 监测 加 载 过 程 中 的 应 力 应 变 
以 及 泊 松 比 的 变化 ,提取 最 大 单 轴 抗 压强 度 o, ,并 
以 o./2 处 的 切线 模 量 作为 弹性 模 量 ; 同 理 ,巴西 
辟 裂 数值 试验 以 应 变 率 0.5 x 107/s 控制 进行 试 
验 。 而 在 进行 含 结构 面试 件 的 剪 切 试验 时 , 先 对 试 
样 施加 正 应 力 ,然后 设 定 对 剪 切 盒 施加 剪 切 应 变 率 
为 0.2。 
以 室内 岩石 力学 试验 参数 和 试 件 破 坏 类 型 为 依 
据 ,通过 不 断 调整 细 观 参数 进行 模拟 ,最 终 获 得 了 与 
室内 试验 相 匹 配 的 宏观 力学 参数 ( 表 1)、 应 力 应 变 
曲线 (图 5) 和 试 件 破 坏 类 型 (图 6)。 
20[ -。 室内 试验 


一 一 数值 试验 


0=173.73 MPa 


An 
be 


E=21.25 GPa 


0=171.10 MPa 


We 
SS 


een 


轴 向 应 力 c./MPa 
站 一 


en 
[= 


E=21.55 GPa 


0 0.005 0.010 0.015 
轴 向 应 变 e 


图 5 岩石 数值 匹配 试验 应 力 应 变 曲线 图 
Fig.5 Stress strain curve of numerical matching test 
表 1 室内 试验 与 数值 试验 宏观 力学 参数 对 比 
Tab. 1 


Comparison of macro mechanical parameters 


between laboratory and numerical test 


力学 参数 及 单位 室内 试验 数值 试验 
抗 压 强度 o./MPa 171.1 173.74 
弹性 模 量 已 /GPa 21.55 21. 25 
泊 松 比 v 0.19 0.19 
抗 拉 强度 o,/MPa 11.30 11.31 
强度 比 eever 15.14 15.36 
内 聚 力 C /MPa 0.45 0.44 
内 摩擦 角 pg /(°) 39.00 39. 00 


由 此 可 知 ,以 细 观 力学 参数 ( 表 2) 得 到 的 数值 
试验 宏观 力学 参数 与 岩石 试 件 室内 试验 与 相近 ( 参 
数 误差 均 在 上 0. 022 以 内 ) ; 单 轴 压 缩 条 件 下 的 应 力 
应 变 曲 线 图 中 ,数值 试验 中 岩石 为 理想 弹 脆性 ,而 实 
际 试验 中 压 密 阶段 的 单 轴 压 缩 和 巴西 辟 裂 试验 均 呈 
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现 相 同 的 脆性 拉 破 坏 形式 , 直 剪 试验 试 件 呈 相同 的 ”因此 ,颗粒 流 细 观 力学 参数 ( 表 2) 能 够 反映 深 蚀 兰 
剪 切 破坏 形式 ,表明 模拟 结果 均 与 室内 试验 一 致 。 ” 体 的 力学 行为 。 


[ 颜 


色 ] 不 同 的 独立 破坏 块 体 


脆性 拉 破 坏 剪 切 破坏 
(a) 室内 试验 (b) 数值 试验 
图 6 室内 试验 与 数值 破坏 形式 对 比 
Fig.6 The comparison between the laboratory test and the numerical test 
表 2 岩 体 颗粒 及 其 接触 模型 细 观 参数 


Tab.2 Mesoscopic parameters of rock mass particles and contact model 


类 型 参数 单位 接触 范围 值 
颗粒 半径 比 Ra/ Rin 四 全 部 颗粒 1.66 
最 小 半径 Ri cm 0. 66 
摩擦 系数 = 0. 60 
改进 PBM 杨 氏 模 量 玉 GPa 全 部 颗粒 25.0 
刚度 比 K 2.0 
软化 系数 B 14.0 
内 摩擦 角 gp (°) 58.0 
抗 拉 强 度 o MPa 方解石 6.72 
石英 7.90 
黏土 矿物 4.74 
方解石 -石英 6.32 
方解石 -黏土 4.35 
石英 -黏土 3.95 
内 聚 力 C MPa 方解石 93.5 
石英 110.0 
黏土 矿物 66.0 
方解石 -石英 88.0 
方解石 - 黏 士 矿物 60.5 
石英 -黏土 矿物 55.0 
SJM 法 向 刚度 K， GPa. m-! 结构 面 附近 30.0 
切 向 刚度 K、 GPa .mr-l 20.0 
摩擦 系数 = 0.6 
内 摩擦 角 wm C35 45.0 
剪 胀 角 更 (9 5.0 
抗 拉 强度 or MPa 0.1 
内 聚 力 C， MPa 0.15 
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3 ” 溶 蚀 岩 体 单 轴 压 缩 试验 模拟 


3.1 岩 体 结构 面 参数 统计 


研究 区 内 的 岩 体 结构 面 通常 为 1 组 层面 和 2 组 
节理 面 。 通 过 布设 测 线 和 测 窗 法 对 研究 区 溶 蚀 岩 体 
结构 面 特征 参数 进行 测量 和 统计 。 对 约 300 条 结构 
面 采 用 概率 统计 法 分 析 发 现 ,该 结构 面 倾向 和 倾角 
主要 服从 正 态 分 布 ,而 迹 长 和 间距 服从 负 指 数 分 布 
( 表 3)。 

表 3 溶 蚀 岩 体 结构 面 参 数 


Tab.3 Discontinuity parameters of dissolution rock mass 


正 态 分 布 负 指 数 分 布 
结构 面 
倾向 /(") ”倾角 /(°) 迹 长 /m ”间距 /m 
层 理 面 215 +0 15 +0 - 0.16 +0.08 


节理 1 189.90 +18.97 84.26 +2.83 0.16+0.10 0.15 +008 


节理 2 305.75 +15.92 86.46 +3.46 0.16+ 上 0.06 0.22 +0.02 


3.2 等 效 溶 蚀 岩 体 模型 建立 
由 于 岩 体 具有 斥 才 效 应 , 需 确 定 其 表征 单元 体 


积 (representative elemental volume, REV ) 。 据 结构 
面 参数 及 分 布 ( 表 3) , 借助 裂 际 建 模 软件 将 不 同 尺 
才 的 结构 面 网 络 耦 合 于 颗粒 流 中 ,建立 等 效 节理 岩 
体 进 行 力 学 试验 模拟 。 当 岩 体 力学 性 质 趋 于 稳定 
时 , 即 可 确定 该 最 小 表征 单元 体 的 尺寸 。 数 值 试验 
表明 ,研究 区 的 REV 约 为 1mx1mx2m( 图 7a)。 
基于 岩 体 表征 单元 体 ,采用 与 匹配 试验 同样 的 方法 
和 参数 建立 颗粒 模型 (图 7b) ,并 赋予 各 接触 模型 ， 
建立 等 效 节理 岩 体 模型 (图 7c) 。 

描述 岩溶 ( dissolution ) 发 育 程度 可 以 采用 体 溶 
蚀 率 表示 。 本 研究 溶 刨 率 让 为 体 溶 刨 率 ,为 溶 蚀 
空洞 所 占 体 积 (W) 与 总 体积 (V) 之 比 


V, 
k = x 100% (10) 


溶 蚀 孔洞 采用 自动 元 胞 机 方法 *| 生 成 , 即 在 模 
型 中 的 方解石 颗粒 中 投放 一 定 比例 的 随机 种 子 , 再 
通过 聚 类 算法 在 种 子 颗 粒 周 于 随机 捕获 颗粒 ,最 后 
将 这 些 颗 粒 进行 溶 蚀 剔除 ,分 别 建立 深 蚀 率 上 为 
0% 5% 、10% 15% .20% 和 25% 的 6 种 等 效 的 孔隙 
型 溶 蚀 岩 体 模 型 ,如 图 7(d) 所 示 。 
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Ga) 结构 面 网 络 模型 


i 
0 RE 


(0) 节理 岩 体 接触 模型 


(d) 溶 蚀 岩 体 颗粒 模型 
司 7 等 效 洲 仔 岩 体 模型 


Fig.7 Equivalent model of dissolution rock mass 


3.3 单 轴 压缩 数值 试验 


首先 ,在 岩 体 模型 顶部 和 底部 创建 加 载 板 ,加载 
板 与 溶 蚀 宕 体 模型 之 间 的 摩擦 系数 为 0, 以 消除 试 
验 的 端面 摩擦 ;上 下 加 载 板 应 变 率 为 0.1; 其 次 ,在 
考虑 试 件 端 部 局 部 应 力 集中 效应 后 ,利用 fish 实时 
监测 和 记录 试 件 应 力 应 变 等。 其 中 , 试 件 的 轴 向 应 
力 为 所 受 的 总 接触 力 除 以 面积 , 轴 向 应 变 为 轴 向 位 
移 与 试 件 初始 高 度 之 比 ;最 后 , 当 轴 问 应 力 降低 至 峰 
值 应力 80% 时 ,停止 加 载 。 


4 试验 结果 与 讨论 


4.1 溶 蚀 岩 体 细 观 微 裂 纹 演 化 特征 


岩 体 力学 强度 是 通过 颗粒 间接 触 模型 反映 , 即 
由 平行 黏 结 模型 和 光滑 节理 模型 来 分 别 反映 颗粒 间 
的 胶结 , 岩 体 内 平行 黏 结 的 断裂 可 认为 是 岩石 内 微 
裂纹 的 产生 ,而 节理 符 结 的 断裂 可 认为 是 结构 面 微 
裂纹 的 产生 与 拓展 。 

在 单 轴 压 缩 条 件 下 ,不 同 浴 蚀 率 ( 丰 ) 岩 体 REV 
的 岩 块 和 结构 面 累计 微 裂 纹 数目 随 轴 向 应 变 s 的 
变化 关系 如 图 8 所 示 。 由 图 可 知 ,在 单 轴 加 载 前 期 ， 
结构 面 微 裂 纹 产 生 或 扩展 数量 急剧 增加 ,而 岩 块 内 
部 的 微 裂纹 缓慢 增加 ; 随 着 轴 向 应 变 不 断 增 大 , 岩 块 
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内 部 微 裂 纹 呈 指数 迅速 增加 ,而 结构 面 趋 于 稳定 ,这 
表明 在 单 轴 压 缩 条 件 下 , 岩 体 内 结构 面 首 先 破坏 , 随 
后 转 为 岩 块 内 部 微 裂 纹 损伤 产生 ,进而 导致 涯 块 内 
部 强度 降低 。 

由 图 8(a) 可 知 ,相同 溶 蚀 率 时 ,省 块 内 部 的 
拉 裂 纹 数 ,的 数量 级 为 10 ,是 前 切 裂 纹 数 NN, 的 10 
~20 倍 , 说 明 岩 块 以 脆性 拉 破 坏 为 主 ,损伤 体现 为 
拉 和 裂纹 产生 和 拓展 ;相同 轴 向 应 变 时 , 随 着 溶 蚀 率 的 
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增加 , 拉 裂 纹 数量 越 大 , 增 速 越 快 ,而 剪 切 微 裂纹 变 
化 无 明显 规律 ,表明 溶 蚀 主要 影响 宕 块 内 部 的 拉 裂 
纹 产 生 与 扩展 。 由 图 8(b) 可 知 ,相同 溶 蚀 率 上 时 ， 
剪 切 微 裂 纹 数 NN 为 10 数量 级 ,是 拉 和 裂纹 数 N, 约 50 
信 , 表 明 结 构 面 破坏 以 剪 破 坏 为 主 ,产生 大 量 剪 切 微 
裂纹 ;相同 的 轴 问 应 变 条 件 下 , 剪 切 微 裂纹 随 溶 蚀 率 
的 增加 而 减少 ,而 拉 破 坏 则 也 表现 出 无 明显 规律 , 表 
明 宕 体 结构 面 的 剪 切 裂纹 形成 受 溶 蚀 作用 的 影响 。 
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(b) 结构 面 内 微 裂纹 数 


到 8 溶 蚀 岩 体 中 岩 块 和 结构 


i 微 橡 纹 数 变 化 


Fig.8 Variation of the number of microcracks of rock blocks and discontinuities in dissolution rock mass 


4.2 溶 蚀 岩 体 宏观 变形 破坏 特征 


单 轴 压 缩 条 件 下 ,不 同 溶 蚀 率 的 孔隙 型 岩 体 应 
力 - 应 变 曲 线 与 溶 刨 率 k = 5% 的 累计 微 裂 纹 总 数 入 
关系 如 图 9 所 示 。 由 图 可 知 , 随 溶 蚀 率 增加 , 岩 体 的 
峰值 强度 和 峰值 应 变 不 断 下 降 。 

岩 体 变形 破坏 由 结构 面 和 岩 块 变形 破坏 体 
现 ”“ , 当 贿 体 结 构 面 并 不 能 直接 形成 贯通 性 破坏 面 
时 , 岩 体 的 软弱 结构 面 黏 结 首先 断裂 而 使 得 微 裂 纹 
激增 , 即 结构 面 快 速 破坏 阶段 (OA) ;随后 岩 块 发 生 
弹性 变形 ,内 征 裂 纹 几 乎 不 变 , 岩 体 应 力 -应 变 呈 一 
定 线性 关系 ,为 岩 块 线 弹 性 变形 阶段 (AB) ;此 后 岩 
块 届 服 发 生 塑性 变形 , 岩 体 进入 岩 块 塑性 变形 阶段 
(BC ) ;最 后 , 当 应 力 超过 单 轴 抗 压强 度 时 , 岩 体 进入 


完全 破坏 阶段 (CD ) 。 


120 本 
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图 9 孔隙 型 洲 蚀 岩 体 的 应 力 -应 变 曲线 
Fig.9 Stress-strain curves ofdissolution rock mass 
岩 体 破坏 以 岩 体 分 裂 为 块 体 和 形成 贯通 破坏 面 
为 标志 " 。 由 单 轴 压 缩 试验 的 溶 蚀 岩 体 的 破坏 形态 
特征 图 (图 10) 可 知 , 碎 块 数量 大 、 尺 寸 小 , 微 裂 纹 遍 
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布 岩 体 ,表明 岩 体 各 部 分 受 力 均匀 ,能 耗 散 较 多 能 
量 ; 随 着 浴 蚀 率 的 增加 ,破坏 碎 块 数量 减少 , 碎 块 尺 
寸 增 大 , 微 裂纹 呈 集 中 分 布 , 岩 体 中 因 单 轴 加 载 产 生 
的 裂纹 集中 于 洞 际 附近 ,导致 竺 体 从 均匀 性 破坏 的 
转化 为 局 部 的 结构 性 (结构 面 和 洲 蚀 洞 陈 等) 破坏 ， 
进而 形成 贯通 性 破坏 面 ,降低 了 岩 体 的 强度 。 
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(b) 微 裂纹 分 布 
图 10 不 同 溶 蚀 率 涯 体 的 破坏 形态 


Fig. 10 Failure modes of rock mass with 


different dissolution rates 
4.3 溶 蚀 率 与 岩 体 力学 强度 参数 的 关系 


目前 ,通过 室内 岩石 力学 试验 方法 获得 溶 蚀 涯 
体 强度 参数 存在 很 大 困难 。 因 此 ,通过 建立 溶 蚀 率 
与 岩 体 强度 参数 的 关系 ,是 解决 洲 蚀 岩 体 强度 参数 
获取 困难 的 重要 途径 。 

单 轴 抗 压强 度 和 变形 模 量 是 两 个 常用 的 参数 ， 
通过 单 轴 压 缩 条 件 下 不 同 溶 蚀 率 的 力学 参数 试验 ， 
获取 单 轴 抗 压强 度 o, 与 岩 体 变形 模 量 5, ,最 后 拟 合 
得 到 溶 蚀 率 与 抗 压强 度 和 变形 模 量 的 函数 关系 (图 
11)。 由 图 可 知 , 溶 蚀 发 育 程度 较 低 时 ,引起 岩 体 的 
单 轴 抗 压强 度 降低 幅度 较 大 ; 当 涂 人 蚀 率 不 断 增加 ,其 
降低 幅度 不 断 减 小 , 该 规律 与 朱 雷 等 ” 利用 
FLAC ”研究 溶 蚀 砾 岩 得 到 的 抗 剪 强度 与 溶 蚀 程度 关 
系 基 本 一 致 。 孔 隙 型 浴 蚀 岩 体 的 单 轴 抗 压强 度 or 
和 游人 蚀 率 丰 呈 负 指数 关系 , 即 


ae = aoe (11) 
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式 中 ,au 为 节理 岩 体 单 轴 抗 压强 度 ,单位 为 MPa。 

随 着 浴 蚀 发 育 程度 不 断 增 大 , 溶 蚀 岩 体 的 变形 
模 量 呈 线 性 降低 ,斜率 为 -17.54, 即 变形 模 量 ,与 
洲 蚀 率 正 呈 反 比例 关系 为 


E, = Eo -17.54k (12) 
式 中 ,io 为 节理 岩 体 变 形 模 量 , 单 位 为 GPa。 
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11 溶 蚀 率 与 溶 蚀 岩 体 力学 参数 的 函数 关系 
Fig.11 Function relationship between dissolution rate and 


mechanical parameters of the rock mass 


5 结 论 


本 研究 以 贵州 喀斯特 地 区 溶 蚀 岩 体 为 研究 对 
象 ,基于 岩石 矿物 含量 和 代表 性 岩 体 结构 面 特征 , 利 
用 三 维 颗粒 流离 散 元 法 研究 单 轴 压 缩 条 件 下 孔隙 型 
溶 蚀 岩 体 的 力学 特性 ,得 出 以 下 结论 。 

1) 采 用 改进 的 平行 忒 结 模 型 ,基于 宕 石 细 观 矿 
物 特征 和 表征 单元 体 结构 面 网 络 ,并 利用 多 宏观 参 
数 标定 可 建立 等 效 节 理 岩 体 模型 ;在 节理 岩 体 模型 
中 采用 元 胞 自动 机 算法 将 部 分 溶 蚀 颗粒 剔除 ,可 建 
立 随机 溶 蚀 岩 体 模型 。 

2) 单 轴 加 载 初期 , 岩 体 内 结构 面 剪 切 微 裂纹 数 
激增 ,而 岩 块 内 微 裂纹 为 零 ; 随 着 应 变 逐 渐 增 大 , 岩 
块 内 以 拉 裂 纹 为 主 , 微 裂纹 数 呈 指数 增加 。 

3 ) 随 着 洲 蚀 率 增 加 , 岩 体 内 微 裂纹 总 数 不 断 降 
低 ; 相 同 轴 向 应 变 时 , 随 溶 蚀 率 增加 , 岩 块 内 拉 伸 微 
裂纹 增加 , 而 结构 面 剪 切 微 裂纹 减少 。 

4) 岩 体 应 力 -应 变 曲线 可 分 为 结构 面 快 速 破 坏 、 
岩 块 弹性 变形 . 岩 块 塑性 变形 和 岩 体 完全 破坏 等 4 
个 阶段 ; 随 涂 蚀 率 增加 ,峰值 强度 和 和 峰值 应 变 不 断 
降低 , 且 岩 体 从 均匀 性 破坏 转 为 局 部 结构 性 破坏 。 

5) 孔 隙 型 溶 蚀 岩 体 单 轴 抗 压 强度 和 变形 模 量 参 
数 与 溶 蚀 率 分 别 呈 负 指 数 和 反比 函数 关系 ,可 用 于 
估算 溶 蚀 岩 体 的 力学 参数 ,这 为 解决 游人 蚀 岩 体 强度 
室内 试验 难以 获取 的 问题 提供 新 的 途径 。 
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